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Abstract. The present constraints on fundamental constant variation (a and /i) obtained in 
the radio range are reviewed, coming essentially from absorption lines in front of quasars 
of intermediate to high-z galaxies, through CO, HI, OH, HCO + , HCN .. lines up to NH 3 
and CII. With ALMA, the sensitivity to detect radio continuum sources in narrow bands 
will increase by an order of magnitude, and the expected progress is quantified. The relative 
advantage of the radio domain with respect to the optical one is emphasized. 
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1. Introduction 

Molecular absorption lines in the radio domain 
are very narrow, and the high spectral resolu- 
tion is well suited to the study of variation of 
constants, with space and time. Observations 
of absorption lines in the millimeter domain 
have been developped in the 1990's in both 
our own Milky Way (Lucas & Liszt 1994) and 
in very distant galaxies (Wiklind & Combes 
1994). Up to now, at intermediate and high 
redshift, between z=0.25 to z=0.89, only four 
absorption lines systems have been detected 
in the millimeter range and a fifth system at 
0.765, at the OH- 18cm lines (Kanekar et al 
2005). Out of these 5 systems, 3 are inter- 
vening lensing galaxies (and the background 
quasar is multiply imaged), and 2 correspond 
to an absorption of the host (PKS1413+135, 
B3-1504+377). Provided that the background 
continuum source is strong enough, a large 
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number of different molecules and rotational 
transitions can been observed (for an overview 
see Combes & Wiklind 1996; Wiklind & 
Combes 1994, 1995, 1996a,b, 1997). The de- 
tectability depends mainly on the background 
source intensity, and with the increased sensi- 
tivity of ALMA, we could expect to detect 30- 
100 times more sources. 



2. Space variation at z=0 

Even at z=0, molecular emission and absorp- 
tion lines bring very stringent constraints on 
the variation of constants with space. In the 
current quintessence models, proposed to solve 
the dark energy problem, hypothetical scalar 
fields couple to the ordinary matter, leading 
to long-range forces and variations of the fun- 
damental constants (e.g. Garcia-Berro et al 
2007). In particular, the Chameleon theory 
predicts different coupling constants, and dif- 
ferent baryon masses between low and high 
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densities, or for different gravitational poten- 
tials. Since the difference in densities between 
the Earth environment, and molecular clouds 
can be as large as 20 orders of magnitude, 
the Milky Way itself is an ideal laboratory 
to test such variations. One of the abundant 
molecules, which has the largest sensitivity to 
the electron to proton mass ratio fi = m e /m p , 
is the ammoniac, NH3, and its peculiar inver- 
sion lines (Flambaum & Kozlov 2007). It can 
be shown that the position of the lines are 
Av/c - Ki A/i//i, with A", = 1 for rotation lines, 
and Ki- 4.46 for NH3, therefore comparing the 
two give (zinv - z ro ,)/(l + z) = 3.46 Ayu/jU. 

Compared to rotational lines of CCS, 
HC3N, with different sensitivity to yu, it was 
possible to constrain the variation Afi//i at 
the 10~ 8 level. There might remain system- 
atic kinematical variations, differences of AV 
= 10 m/s corresponding to A////i = 10~ 8 
(Levshakov et al 2008, Molaro et al 2009). 

3. High redshift absorption 

Then to probe time variations, a large number 
of high redshift systems are required. However, 
these are quite rare in molecular absorptions, 
since dense molecular clouds are found only 
in the disk of spiral galaxies, and the cross- 
section is small in front of background quasars. 
The first system was found towards the BL Lac 
object PKS1413+135 (Wilkind & Combes 94). 
It is a self-absorbing system, where HI-21cm 
absorption was already detected, pointing to 
high column densities. When an intervening 
galaxy is on the line of sight towards a quasar, 
the probability of lensing is high, and the am- 
plification of the quasar flux is a selection bias. 
This is therefore not a surprise that three out 
of the five known systems are gravitationally 
lensed objects: PKS 1830-21 1, B021 8+357 and 
PMN JO 134-0931. Due to the special geome- 
try for optimum amplification, intervening red- 
shifts range up to z~l, while the quasar is at 
z~2. It becomes difficult to find higher red- 
shift quasars, that are strong enough in the mil- 
limeter, since most radio sources have a syn- 
chrotron steep spectrum. The cosmological K- 
correction works against us in this domain. A 
solution is to follow the higher redshift contin- 



uum sources to the centimetric domain. Let us 
briefly describe the specific features of the ab- 
sorbing systems. 

3.1. PKS1413+135z=0.247 

This is the narrowest absorbing system, with 
two velocity components, each narrower than 
1 km/s. It is a BL Lac, very variable contin- 
uum source, since the radio jets are seen al- 
most along the line of sight. This small angle, 
with a large Lorentz factor of the jets, produces 
superluminal motions. The absorbing galaxy 
appears edge-on, and the absorbing molecular 
cloud is CO optically thick, with N(H 2 ) > 10 22 
cirT 2 , and an visual extinction of A,, > 30 mag. 
The comparison of the position of the HI-21cm 
line and molecular lines constrains the varia- 
tion of constants (Wiklind & Combes 1997). 

3.2. B3 1504+377 z=0.672 

This is the second self-absorbing system, with 
a narrow component (may be corresponding 
to the outer disk), and a broad component, at 
330km/s (may be widened by non-circular mo- 
tions in the nucleus). Up to 7 different molecu- 
lar lines have been detected, and the chemistry 
of the absorbing diffuse gas is unexpected; with 
high HCO + abundance. The comparison of the 
HI-21cm position with the molecular lines re- 
veals strong differences in velocity, which pre- 
vent to use this system for the constant varia- 
tion study. 

3.3. B021 8+357 z=0.685 

This system corresponds to the largest col- 
umn density N(H 2 ) ~ 10 24 cm" 2 ; all three CO 
isotopes (CO, 13 CO and C ls O) have optically 
thick lines, and this allowed the detection of 
many molecules (Combes & Wiklind 1997). It 
was for instance possible to search for molecu- 
lar oxygen O2 without atmospheric absorption, 
with the redshifted lines at 368 and 424 GHz: 
the abundance ratio was constrained to O2/CO 
< 2 10~ 3 , meaning that most of the oxygen 
should be in atomic form (Combes & Wiklind 
1995, Combes et al 1997). 
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The system is a gravitational lens, with two 
images of the quasar, and an Einstein ring. The 
image separation is 335mas (1.8kpc), well de- 
termined through radio interferometry (VLA, 
MERLIN, Browne et al 1993). Substructures 
have been resolved in each of the two images, 
in the form of an extended radio jet compo- 
nent and two bright cores, through VLBI at 8.4 
GHz, with lmas synthetized beam (Biggs et al 
2003). The lens is a spiral galaxy almost face- 
on, and spiral arms complicate the lens analy- 
sis, to obtain the Hubble constant value from 
the time delay (York et al 2005). 

3.4. PKS1 830-211 z=0.88582 

This system is also a gravitational lens with 
two images and an Einstein ring. The redshift 
of the lens was discovered in the millimeter 
range, through sweeping of the band with the 
SEST telescope (after 14 GHz sweep in 14 tun- 
ings, already 2 lines were detected, Wiklind & 
Combes 1996a, 1998). The two images were 
resolved with OVRO by Frye et al (1997). 
There are two velocity components, each one 
absorbing in front a different quasar image. 
This allows to monitor the variations and the 
time delay between the two images with a sin- 
gle dish only, without resolving the 2 images. 
The monitoring during 3 years (lh per week) 
determined a delay of 24+5 days, yielding a 
value for the Hubble constant of Ho = 69 + 12 
km/s/Mpc. An isotopic survey of the various 
molecules revealed different abundances than 
at z=0, due to the predominance of young stars, 
given the relative youth of the galaxy (Muller 
etal2006). 

3.5. PMNJ0134-0931 z=0.7645 

The z=0.7645 lens produces 6 images of the 
background quasar at z=2.22. HI-21cm and 
OH- 18cm lines have been detected at GBT 
(Kanekar et al 2005). The main OH lines at 
1665 and 1667 MHz are in absorption, while 
the satellite line at 1720 MHz is in emission. 
The other satellite line at 1612 MHz is in ab- 
sorption, such that the sum of the two satellite 
lines is compatible with noise. The two satel- 



lite lines are therefore conjugate, which comes 
from the peculiar excitation of the upper levels 
through far-infrared radiation. This conjugate 
character ensures that the gas comes from the 
same regions, and therefore intrinsic velocity 
offsets will not perturb the constant variation 
study. 

4. Variations of constants 

Constraining results have been obtained first 
in the optical domain, through absorption lines 
in front of quasars, with the Alkali Doublet 
(AD) method (lines of CIV, Sill, SilV, Mgll, 
A1III...), and then the "many multiplet" method 
which takes advantage of a greater statistics 
(see e.g. the review by Petitjean et al 2009). 
A possible variation of the fine structure con- 
stant Aa/a = (-0.6 ± 0.15) 10~ 5 was reported 
(Murphy et al 2003) over redshifts between 
0.2 and 3.7. However, with the same method 
at VLT, Srianand et al (2004) and Chand et al 
(2006) do not find any variation, with Aa/a = 
(-0.6 + 0.6) 10~ 6 over a comparable redshift 
range. It is therefore interesting to compare 
with other methods. 

4.1. Constant variation with radio 
absorption 

According to the lines and molecules consid- 
ered, it is possible to have constraints on a, 
fi = m e /m p , the gyromagnetic constant of the 
proton g p , or a combination of them, such as 
y = a 2 g p fi. The HI 21cm hyperfine transi- 
tion is proportional to n p /j. B /ha 3 , with the Bohr 
radius a = h 2 /m e e 2 , and \x p = g p eh/4-m p c, 
yUg = eh/2m e c. On their side, the CO rota- 
tion lines frequencies vary as h/ (Ma 2 ), with M 
the reduced mass of the molecule. The great 
advantages of the radio domain are the large 
spectral resolution, up to R = 10 6 thanks to 
the heterodyne techniques, and the very narrow 
lines due to cold gas on the line of sight. 

4.2. Variation of y, a and fj. 

A global comparison of all molecular lines ob- 
served with the HI-21cm absorption lines in 
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PKS1413 and B0218 systems, the two nar- 
rowest line systems, have given quite strin- 
gent constraints, Ay/y = (-0.20 ± 0.44) 10" 5 
for PKS1413 and Ay/y = (-0.16 + 0.54) 10" 5 
for B0218 (Murphy et al 2001). The preci- 
sion is comparable to the MM method, with a 
limited number of absorbing systems. The ra- 
dio method would then be very powerful with 
more systems detected. The main problem with 
this method is a possible kinematical bias, i. 
e. that the HI absorbing medium is not at the 
same velocity as the CO absorbing one. Of 
course, it is an issue only if a positive vari- 
ation is detected. Since the velocity offset is 
expected to be randomly distributed, the so- 
lution is to observe a large number of sys- 
tems. The precision will then be only limited 
by the dispersion in the offsets. This can be 
measured on the large number of HI and HCO + 
absorbing lines detected towards quasars in the 
Milky Way (Lucas & Liszt 1998). The dis- 
persion is of only 1.2 km/s, corresponding to 
Ay/y = 0.4 10~ 5 . 

In the special system PMN0134, although 
no HCO + nor other molecules have been de- 
tected with IRAM, ATCA or VLA, the use 
of the OH lines, together with HI, have also 
brought interesting constraints, less dependent 
on a possible velocity offset, because of the 
conjugate emission/absorption of the satellite 
OH lines. Together with B0218, the value of 
F = g p [a 2 /jj] 151 has been found not to vary, 
i.e. AF/F = (0.44 + 036 stat + 1.0^0 10 -5 
for < z < 0.7 (Kanekar et al 2005). This 
is a null variation with a sensitivity at 2<x of 
Aa/a ~ 6.7 10" 6 and Afi/ji ~ 1.4 10~ 5 over 
half of the age of the universe. 

4.3. Radio 21cm and UV combined 

Another method has combined several systems 
observed both in HI-21cm absorption and in 
UV, to constrain the parameter y = a 2 g p ii. 
With 9 absorbing systems, over a redshift range 
0.23 < z < 2.35, null evolution was found, 
with Ay/y = (0.63 + 1) 10~ 5 or dy/dt/y = 
(0.6 + 1.2) 10" 15 /yr (Tzanavaris et al 2007). 
The precision is clearly limited by intrinsic line 
of sight velocity differences, with random sign, 
of amplitude ~ 6 km/s. Combining with a vari- 



ation constraints, and restricting at high red- 
shift systems, this gives a constraint in // of 
Afi/fi = (0.58 + 1.95) 10~ 5 , which is in con- 
tradiction with the results from Reinhold et al 
(2006), who found with Lyman band of H 2 
lines a positive detection of variation, /j. being 
higher in the past. 

4.4. Radio centimeter lines 

The high sensitivity if the NH3 inversion lines 
to variation in the fi ratio was used by Henkel 
et al (2009) in a recent multi-line study of 
PKS1830-211 atz~0.9. 

They have observed with the Effelsberg 
telescope NH 3 , HC 3 N, SiO, C 34 S, H 13 CO + , 
H 13 CN, i.e. isotopic lines which are optically 
thin. While systematic errors are AV < 1 km/s, 
they find a limit of A/j./fi < 1.4 10~ 6 . 

The NH3 radio inversion lines were also 
used by Murphy et al (2008) in B021 8, at z=0.7 
who find no variations (A/i/jU < 1.8 10~ 6 , cor- 
responding to Av < 1.9km/s), and the same 
negative result was found with quasar Lyman 
band H2 absorption lines, with improved spec- 
tral calibration and fitting all lines in the Lya 
forest, by King et al (2008): Ajj/fj. = (2.6 ± 
3.0) 10" 6 with redshift between 2 and 3. 

While fitting the B0218 absorption lines, 
Murphy et al (2008) split the profiles in a 
large number of gaussian components, which 
could sometimes appear unrealistic. However, 
the optimal number of components minimizing 
the^ 2 is 8. 

5. Main caveats with the radio 
molecular absorption method 

Up to now, fortunately, the radio absorbers 
have provided only upper limits, and have not 
too much suffered from these caveats: 

1- The assumption that all absorbing com- 
ponents are at the same velocity, whatever the 
molecule could be wrong, and this might hap- 
pen for several reasons, chemistry, excitation, 
temperature, density etc... 

2- The lines could be at different optical 
depths, which might shift the resulting average 
velocity 
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3- Some lines have resolved hyperfine 
components (HCN, NH3, ..) and the assump- 
tion of LTE for their excitation might fail 

4- There could be completely different ma- 
terial in projection on the line of sight, with dif- 
ferent origin and abundances, or different cov- 
ering factor 

5- At the various frequencies observed, 
the background continuum source has differ- 
ent sizes, and the absorbing medium is not the 
same 

6- The continuum source varies in both in- 
tensity and shape; and the compared lines have 
not been observed simultaneously 

Some of these problems could be cured 
by selecting high density tracers, optically thin 
isotopes, lines at similar continuum frequen- 
cies, and taking care of the simultaneity of ob- 
servations. In any case, for the systematic ve- 
locity offsets, or the different chemistry com- 
ponents, only large statistics will smooth the 
errors down. 

6. Advantages of the radio method 

The radio method has the potentiality to give 
the best constraints, provided that more sys- 
tems are detected, since: 

(1) The rest frequencies are more precisely 
known than in the optical, and there is no prob- 
lem of frequency calibration, at the precision 
level required 

(2) The sensitivity of the line position to the 
variation of constants is much higher (K, larger 
by a factor 100!) 

(3) The isotopic lines are observed sepa- 
rated, while in the optical there is a blend, with 
possible differential saturation 

(4) Also the hyperfine structure is resolved, 
while in the optical there could be differential 
saturation, and velocity shifts 

(5) The spectral resolution can be very high 
(10 6 , heterodyne) and the observed lines are 
very narrow 

7. Progress with ALMA 

The radio studies suffer from the small num- 
bers of sources, and the fact that they are all 
at low redshift. Fortunately, the sensitivity of 



ALMA will be able to detect many more con- 
tinuum sources, to search for absorption lines, 
and at larger redshifts. ALMA will have a 
much wider bandwidth, allowing the search of 
absorption, even if not previously detected in 
the optical or HI-21cm. The redshift will be ob- 
tained directly in the millimeter. 

Up to now, there is about 50 systems de- 
tected in absorption in HI-21cm, at z > 0.1. 
Curran et al (2006, 2007, 2008) have under- 
taken a new search in front of 13 sources, and 
found only one new HI-21cm detection. With 
GMRT, Gupta et al (2009) have found inter- 
vening HI-21cm absorbers, with a higher de- 
tection rate: 9 new 21cm detections out of 35 
Mgll systems, at z~ 1 .3. It is possible to predict 
the number of continuum sources that can be 
selected as targets for absorption searches with 
ALMA. The density of flat-spectrum quasars 
has been shown to follow the same curve as 
optical quasars, a curve peaking at z ~ 2, very 
similar to the star formation history (Wall et 
al 2005). They are however still of significant 
density atz ~ 3. 



8. Conclusion and perspectives 

Due to caveats (velocity offsets, variability, op- 
tical depths, hyperfine ratios, etc) a negative re- 
sult is more easy to believe than a positive one! 

There is many ways to improve the results: 
first obtaining more systems and more statis- 
tics, then improving the sensitivity and the 
signal-to-noise ratios, taking into account the 
simultaneity of observations for the lines to be 
compared. The number of continuum sources 
that will be use as targets for absorption stud- 
ies will be one or two orders of magnitude 
larger than now with ALMA, although the mil- 
limeter domain is not favoured here by the K- 
correction. 

It could then be interesting to search 
rest-frame 3mm systems at centimeter wave- 
lengths, with EVLA, SKA precursors and in 
the future with SKA. 
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Fig. 1. Comparison between all observed constraints on variation for a (top) and ji (down). The data in the 
a diagram are from OKLO (Damour & Dyson 1996), MMM (Murphy et al 2003, Srianand et al 2004), Fell 
(Quast et al 2004) and radio (Murphy et al 2001). The data in the yu diagram are from 21-UV (Tzanavaris et 
al 2007), H, (Reinhold et al 2006, King et al 2008), and radio (Murphy et al 2008, Henkel et al 2009). 



References 



Biggs, A. D., Wucknitz, O., Porcas, R. W. et al.: 
2003 MNRAS 338, 599 



Browne, I. W. A., Patnaik, A. R., Walsh, D., 
Wilkinson, P. N.: 1993 MNRAS 263, L32 

Chand H., Srianand R., Petitjean P. et al.: 2006, 
A&A451,45 

Combes F, Wiklind T: 1995, A&A 303, L61 



288 



Combes: Future of radio measurements 



Combes, R, Wiklind, T. 1996, in Cold Gas at High 
Redshift, eds. Bremer, M., et al. Kluwer, p. 215 
Combes R, Wiklind T.: 1997, ApJ 486, L79 
Combes, R, Wiklind, T., Nakai, N.: 1997 A&A 327, 
L17 

Curran, S. J., Whiting, M. T., Murphy, M. T. et al.: 

2006 MNRAS 371,431 
Curran S., Tzanavaris R, Murphy M. et al : 2007, 

MNRAS 381, L6 
Curran S., Whiting M„ Wiklind T. et al : 2008, 

MNRAS 391, 765 
Flambaum V.V., Kozlov M.G.: 2007, Phys Rev Lett 

98, 240801 

Frye, B., Welch, W. J., Broadhurst, T.: 1997 ApJ 
478, L25 

Garcia-Berro, E., Isern, J., Kubyshin, Y. A.: 2007, 

A&ARevl4, 113 
Gupta N., Srianand, R., Petitjean P. et al: 2009, 

MNRAS 398, 201 
Henkel C, Menten K., Murphy M. et al: 2009 A&A 

500,725 

Kanekar N., Carilli C, Langston G. et al: 2005 

PhRvL 95, z 1301 
King, J. A., Webb, J. K., Murphy, M. T., Carswell, 

R. R: 2008, Phys Rev Lett 101, 1304 
Levshakov S., Molaro P., Kozlov M.: 2008, preprint, 

arXiv:0808.0583 
Lucas R., Liszt, H.S. 1994, A&A 282, L5 
Lucas R., Liszt, H.S. 1998 A&A 337, 46 
Molaro P., Levshakov S., Kozlov M.: 2009, preprint, 

arXiv:0907.1192 
Muller S., Guelin M., Dumle M., Lucas R., Combes 

F: 2006, A&A 458, 417 
Murphy, M. T., Webb, J. K, Flambaum, V. V., 

Drinkwater, M. J., Combes, F, Wiklind, T.: 

2001 MNRAS 327, 1244 
Murphy, M. T., Webb, J. K, Flambaum, V. V.: 2003 

MNRAS 345, 609 
Murphy M., Flambaum V., Muller S., Henkel C: 

2008, Science 320, 1611 
Petitjean P., Srianand R., Chand R. et al.: 2009, 

SSRv. preprint (arXiv:0905.1516) 
Quast, R., Reimers, D., Levshakov, S. A. 2004 A&A 

415, L7 

Reinhold, E., Buning, R., Hollenstein, U., Ivanchik, 

A., Petitjean, P., Ubachs, W.: 2006 Phys Rev 

Lett 96, 1101 
Srianand R., Chand H., Petitjean P., Aracil B.: 2004, 

Phys Rev Lett 92, 121302 
Tzanavaris, P., Murphy, M. T., Webb, J. K, 

Flambaum, V. V., Curran, S. J.: 2007 MNRAS 

374, 634 

Wall, J. V., Jackson, C. A., Shaver, P. A., Hook, I. 

M., Kellermann, K. I.: 2005 A&A 434, 133 
Wiklind, T., Combes, F. 1994, A&A 286, L9 



Wiklind, T., Combes, F. 1995, A&A 299, 382 
Wiklind, T., Combes, F. 1996a, Nature, 379, 139 
Wiklind, T., Combes, F. 1996b, A&A 315, 86 
Wiklind, T., Combes, F. 1997, A&A 328, 48 
Wiklind, T., Combes, F. 1998, ApJ 500, 129 
York, T., Jackson, N., Browne, I. W. A et al: 2005, 
MNRAS 357, 124 



